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Abstract 

The investigations on the flux relaxation phenomenon of a type-II superconductor are important 
because they provide the information about the flux pinning ability and current-carrying ability 
of the superconductor. However, a unified theory of flux relaxation is currently unavailable. Here 
I present a general mathematical model of the flux relaxation. In this model, I proposed a series 
expansion to the activation energy and derived a general formula for the current decay behavior. 
In the light of these formulas, I can analyze the experimental data on the current decay behavior 
and then calculate the activation energy of a vortex system without subjecting to any special 
conditions. The results are accurate for the current decay measurements from a -B^.S^CaC^Os+x 
superconductor. 

PACS numbers: 74.25.Wx, 74.25.Uv 
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I. INTRODUCTION 



The persistent current (or trapped magnetic field) in a type-II superconductor decays be- 
cause the vortex lines jump between adjacent pinning centers spontaneously due to thermal 
activation^, quantum tunneling^ or mechanical vibration. Flux relaxation has dependency 
on activation energy and temperature. The activation energy is a function of current because 
the Lorentz force of a current reduces the activation energy. By proposing a detailed current 
dependent activation energy, one can obtain the corresponding current decay model. In 
low temperature superconductors, the strong pinning dominates and current decay usually 
shows logarithmic behavior.- This can be explained by Anderson's flux creep theory in which 
he assumed an activation energy with linear current dependence. 6 However, in high temper- 
ature superconductors, the weak point pinning dominates and current decay usually shows 
non-logarithmic behavior. 5 To address this phenomenon, a number of theoretical models 
were proposed, such as collective creep theory^ and vortex glass theory^, both of which give 
an activation energy with nonlinear current dependence*^. Other theoretical models were 
also proposed for the non- logarithmic behavior of current decay, like bulk-surface pinning 
transition^ - — and flux redistribution^. 

The work described here is motivated by the fact that the flux relaxation phenomenon is 
currently described with multiple models. The experimental confirmation of these theories 
are limited to the special cases.— - — In practice, the temperature and current are arbitrary. 
The discrepancy between the theoretical predications and experimental results are large. 
Thus, it is highly desirable to search for a general theoretical model which gives an arbi- 
trary description of the activation energy and corresponding current decay behavior without 
subjecting it to any restraint conditions. The studies on flux relaxation are important be- 
cause the flux relaxation influences the current carrying ability and stability of a type-II 
superconductor; consequently, it limits the possible applications of the superconductor. The 
information about the flux relaxation is also important in understanding the pinning mech- 
anism and vortex structure in a type-II superconductor. 

In this paper, I investigated the general expressions of flux relaxation phenomenon from 
a mathematical framework. First, I expanded the activation energy as a Taylor series of 
current density and obtained the corresponding flux relaxation model using the inverse series 
methods. Next, I discussed the possible physical meaning of various parameters. Finally, 
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I studied the applied conditions of the quadratic activation energy and linear activation 
energy. 

II. MODEL 

Flux relaxation can be described by hopping rate, R cx e~ Ua ' kT , where U a is the activation 
energy of vortices, k is Boltzmann constant and T is temperature.- The hopping rate of the 
vortices is increased by Lorentz force fi = j x B, where j is current density and B is 
magnetic field. Therefore, U a is a decreasing function of j.(Ref. \s\) The thermal hopping 
of vortices causes a fall in j, and j is then a function of time. By proposing a detailed j 
dependent U a (j), one can obtain the time evolution of j(t). This suggests that an important 
step in studying the flux relaxation phenomenon is to search for a detailed j dependent 
U a (j). Based on different physical considerations^"— >2£, physicists obtained a number of 
U a {j). Here I used an approach mainly based on mathematical consideration. 

First, I argue that the activation energy U a (j) must be a nonlinear function of current 
density j because of the elasticity of vortices and interaction between vortices.^S A vortex 
deforms under a driving force. It may yield under a strong driving force and show plastic 
behavior— In this case, the strain energy of a vortex is a nonlinear function of j. On the 
other hand, the interaction energy between vortices is U int = s ® 2 ° x2 K (-0, where r is the 
distance between two vortices, A is the penetration depth and Kq is a zeroth order Hankel 
function with an imaginary argument.- Function K causes a nonlinear response in U a (j) 
with respect to j. Therefore, the U a (j) of a real vortex system must be a nonlinear function 
of j. 

The exact form of the activation energy is usually unknown because of the complexity 
of the interaction between vortices. From mathematics we know that series expansion is a 
good approach for expressing an implicit function. Let us now consider the possibility of 
expanding the activation energy U a (j) as a Taylor series with respect to j. 

For a vortex system, we only need to consider U a (j) on the closed interval [0, j c ]. The 
hopping rate, R oc e~ Ua ^ kT , has limited values because a vortex system decays at any circum- 
stances due to the thermal motion of the vortices. It indicates that U a (j) also has limited 
values; otherwise, a vortex system does not decay. At vanishing driving force, U a (j) is the 
pinning potential U c , i.e., U a (j — > 0) — > U c . U c is a limited number.™ On the other hand, 
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at critical level, U a (j) approaches to zero, i.e., U a (j — > j c ) — » 0. Thus, U a {j) is a bounded 
function for all j G [0, j c ]. Therefore, in most cases one can expand U a (j) as a Taylor series 
of j (at least up to order one, i.e., the linear activation energy). Let us now write out the 
general expression explicitly: 

n 

U a (j) = U c -J2<iif (!) 

i=l 

where U c = U a (0), a x = -U' a {0), a 2 = -E£'(0)/2!, • ■ • , a n = -Ui n) (0) jn\. Normally, the 
activation energy U a has dependence on temperature T and current density j. Further- 
more, U a includes the contribution from the elastic deformation and interaction of vortices, 
which has a dependence on the penetration depth A and coherence length £. Therefore, the 
coefficients in Eq.([T]) should be in the form of U C (T, A, £) and a n (T, A, ^). 

To obtain the current decay model j(t) from Eq.(JTJ), we need the connection between the 
activation energy U a (j) and time t. Early studies^ have shown that, by assuming that flux 
relaxation is caused by thermal activation, the time evolution of U a (j(t)) (with logarithmic 
accuracy) is, 

U a (j(t)) = kTln(l + t/t ) (2) 

where to is a short time scale for the Bean model to be formed.- This is a macroscopic 
quantity depending on sample size. Substituting Eq.© into Eq.flTJ, we have 

n 

w(t) = J2^f(t) (3) 

i=i 

where 

w(t) = U C - kTln(l + t/t ). (4) 

Now our task is to find out the inverse function j(w(t)) from Eq.Q. First, let us expand 
j(w(t)) as a series of w(t), 

n 

j(w(t)) = J2b l w\t) (5) 

i=i 

Employing Cauchy's formula^ from complex analysis, one can easily obtained the 
coefficients b^s 
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k A 1 £ Mr™- (-)'(-)'■■■ 

a%n \aij \aij 

n(n + l)---(n-l + s + t + u-\ ) 

sltlul ■ ■ ■ 

where s + 2t + 3u + ■ ■ ■ = n — 1. From Eq.([6]) we see that only coefficient a\ contributes to the 
logarithmic decay coefficient b±. The coefficients 02,03, ■• ■ ,a n result in the non- vanishing 
coefficients b n (n > 2), which contribute to the deviation away from the logarithmic decay. 

Since Eq.(j3]) and Eq.Q have symmetry, if we switch their sequence, the same procedure 
can be employed to calculate the coefficients a n . Therefore, a n and b n have the same 
symmetry as that of Eq.(|3]) and Eq.(JS]). By simply doing a commutation to the coefficients 
b n «-» a n , we have 

n(n + 1) • • • {n - 1 + s + 1 + u H ) 

sltlul ■ ■ ■ 

where s + 2t + 3m + • • • = n — 1. In practice, we usually first find out the coefficients b n 
from experiments, and then calculate the coefficients a n to obtain U a (j) from Eq.flTJ). The 
U a (j) obtained in this way includes the contributions from the Lorentz force, deformation of 
vortices, interaction between vortices and possibly other unknown sources. Thus, U a (j) is a 
combined response to j. Figure 1 shows that these results are accurate for the experimental 
data of a Bi 2 Sr 2 CaCu 2 8+x single crystal. 

In Eq.(JTJ, we expanded U a (j) as a function of j, the terms including j represent the 
contribution from the driving force. As a result of j being in a superconductor, there will be 
a corresponding magnetic field B according to the Biot-Savart law. So, we can also expand 
U a as a function of the magnetic field B or magnetization M in the same way as we did in 
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Eq.flU, that is, U a {B) = U C -J2 a i B \ or U a {M) = U C -J2 Similar to Eq.©, we can 

i=i i=i 

n 

obtain the time evolution of the magnetic field B(t) = biW l (t), or time evolution of the 

i=l 

n 

magnetization M(t) = ^2biW z (t). 
III. DISCUSSION 

In Eq.©, we expanded the activation energy U a (j) as a Taylor series of current density 
j. Let us now discuss the possible physical meaning of the coefficients. 

1. Meaning of U c . — Eq.© indicates that U a (0) = U c . The constant U c is then the 
activation energy at vanishing driving force, i.e., the pinning potential of a vortex. 

The pinning potential U c is important in determining the irreversibility field Hi = 
_Ec — e u c /k B T w ]j ere e j s the electric field criterion when the current density reduces to 
zero, fiQ is the permeability of vacuum, aj is the flux line spacing and is the attempt fre- 
quency of the vortex.— On the other hand, U c can be used to calculate the critical current 
density j c . From Eq.(T5]) and Eq.(@| we see that the constraint condition U a (j c ) = equals 
to set w(t) = U c . So, we can obtain the critical current density by replacing the w(t) in 
Eq.© with U c , that is, j c = jr biU l c . 

i=l 

2. Meaning of a n (n>2). — From mechanics we know that elastic reaction is represented 
by the quadratic term of independent variables, while inelastic reactions are represented by 
the higher order terms (above two) of independent variables. Let us now prove that the 
activation energy of a elastic vortex system has quadratic current dependence. 

Consider an elastic vortex which is subjected to a uniformly distributed Lorentz force 
fi, the vortex undergoes deformation and adjusts its position to minimize the total energy 
because of the elasticity. For simplicity, let us consider the vortex as a uniform object. The 
elastic strain energy due to the bending of the vortex is U e (j) = j^jdx, where M is the 
internal moment, E is the modulus of elasticity, I is the moment of inertia, L is the length 
of the vortex under consideration and x is the coordinate established along the vortex.— 
Under the uniformly distributed load /l, the equilibrium condition is, M + x (|) /t = 0. 
Assuming that the deformation of the vortex 5 is very small compared with the distance 
between two adjacent pinning centers D (5 << D), one can find the elastic strain energy 
due to the bending of the vortex: U e {j) = = j , where $o is the flux quantum. 
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This suggests that the coefficient a 2 in Eq. ([TJ) represents the weight of the contribution from 
the elastic deformation. 

The above discussion suggests that if we let a n = (n>2) in Eq.flTJ, then it is equal to 
ignore the inelastic deformation and interaction of the vortices. The activation energy of 
the elastic vortices becomes 

U a (j) = U c -a lJ -a 2 j 2 (8) 



substitute Eq.flSJ) into Eq.(J2]), we have (choose one of the solutions which is a decreasing 
function of t), 

1 / 777. 

(9) 



where w(t) is defined by Eq.(j4]). Expanding the right hand side of Eq.([9]) (or using Eq.([6]) 
directly) and keeping the terms up to the second order, we have j(t) ~ -^w(t) — ^w 2 (t). 
The first order term j^wif) represents the logarithmic decay. The second order term ^w 2 (t) 
causes deviation away from the logarithmic decay. Since ^w 2 (t) depends on a 2 , it represents 
the contribution from the elasticity of the vortices. This indicates that any elastic vortex 
under a driving force will show non-logarithmic current decay behavior. 

Replacing the w(t) in Eq.Q with U c , we obtain the critical current density j c = 

3. Meaning of ax. — Further let a 2 = in Eq.® is equal to assume that the vortices 
have very large elastic modulus and elastic deformation can be ignored, i.e., the vortices are 
rigid. The corresponding activation energy becomes 



U a (j) = U C - a x j (10) 

This is Anderson's linear activation energy 6 . Using Eq.(j5]) and Eq.((6]), we obtained the 
logarithmic current model j(t) = —w(t). The corresponding critical current density is 
Jc = -U c . 

Since the Lorentz force on vortices is proportional to the current density j, the linear 
activation energy U a (j) in Eq.f JTOj) is a result of the Lorentz force causes a direct reduction 
of the value a>\j in the pinning potential U c . The early studies have shown that constant % 
is of the order $s<i 2 , where is the total flux in the flux bundle and d is the length of the 
bundled 
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From the above discussion we know that the linear activation energy and corresponding 
logarithmic decay behavior only occurs in the rigid vortex systems. The elasticity and 
interaction of vortices modify the linear activation energy, which cause a deviation away 
from the logarithmic decay behavior. Since the elasticity and interaction of vortices are 
closely related to the penetration depth A of a superconductor,- we can now analyze the 
possible effects of A on the flux relaxation phenomenon. 

The shear modulus Cqq ~ , uniaxial compression modulus and tilting modulus 

Cii ~ C44 ~ |^ 1+ l 2k 2 are all decreasing functions of A, where k is the wave vector.— Low 
temperature superconductors have penetration depth A ~ 500/1, which is smaller than that 
of high temperature superconductors, A ~ (1400 — 2000)A (Ref. |5|) This means that the 
vortices in low temperature superconductors are stiffer, and there should be a better chance 
of observing quasi-logarithmic decays in these superconductors. 

On the other hand, the penetration depth, A(T) « (1 — T/T c ) -1 / 2 , is an increasing func- 
tion of temperature T. The elastic moduli are then decreasing functions of T. At higher 
temperatures, vortices become soft and undergo large deformations, resulting nonlinear ac- 
tivation energy. Therefore, it is more likely to observe non-logarithmic decay behavior at 
higher temperatures, but logarithmic decay behavior at lower temperatures.— 1 ^ 



IV. CONCLUSION 



The activation energy of flux relaxation phenomenon can be expanded as a Taylor series 
of current density. The elasticity of the vortices and interaction between the vortices result 
in the nonlinear activation energy, which causes non-logarithmic decay behavior. It should 
be emphasized that the derivations in this article were mainly based on mathematical skills. 
This ensures a wide application field, but the explicit studies on the physical premises 
should be undertaken according to specific conditions. Secondly, I have neglected the surface 
pinning effects, which may add significant influence on the flux relaxation phenomenon at 
higher temperatures close to the critical level. These need to be considered separately. 
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FIGURE CAPTIONS: 




FIG. 1: (Color online) Time evolution of persistent current. The scattering points are the experi- 
mental data of a Bi 2 Sr2C f aCu%0%+ x (Bi-2212) single crystal, i.e., the persistent current induced at 
27 K under an applied magnetic field of 750 Gauss. 30 The solid black line is the theoretical fit. The 
fitting results are: j(t) = b\w{t) + b 2 w 2 (t) + b 3 w 3 (t) (where w(t) = U c - 27 * k * ln(l + t/t )), b\ = 
(2.77±0.06) x W 5 /k, b 2 = -(1.79±0.01) x 10 3 /k 2 , b 3 = (3.63±0.08)//c 3 , U c = (4.62±0.08) x 10 2 A;, 
t = (1.55 ±0.04). 
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